Revisiting the gravitino dark matter and baryon asymmetry from Q-ball decay 

in gauge mediation 
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We reinvestigate the scenario that the amount of the baryons and the gravitino dark matter is 
naturally explained by the decay of the Q balls in the gauge-mediated SUSY breaking. Equipped by 
the more correct decay rates into gravitinos and baryons recently derived, we find that the scenario 
with the direct production of the gravitino dark matter from the Q-ball decay works naturally. 



I. INTRODUCTION 

The abundances of the baryons and dark matter are 
almost comparable in the present universe. It may be 
understood by the fact that both of them have the same 
origin: Q balls. The Q balls are nontopological soli- 
tons whose stability is guaranteed by the finite charge 
Q [lj]- In supersymmetry (SUSY), they are produced 
during the course of the Affleck-Dine baryogenesis 
so the charge Q is essentially the baryon number. 

In the gauge-mediated SUSY breaking, very large Q 
balls are stable against the decay into nucleons, the light- 
est particle with a unit baryon number. On the other 
hand, if the charge Q is small enough, the Q balls could 
decay into nucleons (baryons). In the latter case, the 
gravitinos are also produced by the Q-ball decay. The 
gravitino is the lightest SUSY particle in the gauge me- 
diation. Therefore, there may be the way to provide both 
the baryons and the dark matter in the universe by the 
Q-ball decay 0413 • 

In Ref. H , two of the present authors investigated the 
scenario to explain both the baryon asymmetry and the 
dark matter in the universe by the Q-ball decay, consider- 
ing the decays into nucleons (quarks), gravitinos, and the 
next to lightest SUSY particles (NLSPs) , taking into ac- 
count the big bang nucleosynthesis (BBN) limits of NLSP 
abundances. Recently in Ref. [IT| . however, another two 
of us studied the Q-ball decay more correctly, and found 
that the decay rate into gravitinos is much smaller than 
that estimated in Ref. [9j. The essence is to consider 
the finite size of the Q ball which becomes smaller than 
the penetration length of the decay particle in the actual 
gauge-mediation Q ball. Therefore, it is necessary to re- 
analyze whether this attractive scenario works or not, 
and this is the purpose of the present article. 

As we will see, it is natural to achieve the right amounts 
of the baryon asymmetry and the gravitino dark mat- 
ter at the same time through the direct production from 
the Q-ball decay although we find the branching into 
gravitinos becomes smaller than the previous estimate 
in Ref. [9]. In this case, more NLSPs are produced by 
the Q-ball decay, but they annihilate afterwards [l(| and 
the remaining NLSP provides only the negligible amount 



of the gravitinos to the dark matter density. Moreover, 
even for m 3 / 2 GeV, severe BBN constraints could be 
avoided if the Q-ball decay temperature is high enough. 

The structure of the article is as follows. In the next 
section, we review the basic features of the Q ball in the 
gauge-mediated SUSY breaking. In Sec. Ill, we provide 
the decay rates of the Q ball, extracted from Ref. [TTj . 
In Sec. IV, we reestimate the abundances of the baryons, 
the gravitinos, and the NLSPs, and seek for the work- 
ing parameters in Scc.V. We conclude in Sec. VI. In the 
appendix, we show, for completeness, some results con- 
cerned with so-called indirect case where the gravitino 
dark matter is produced by the decay of the NLSPs which 
are created by the Q-ball decay. 



II. Q BALL IN THE GAUGE MEDIATION 

The Q ball is constructed from the scalar field called 
the flat direction in SUSY theories. It consists of squarks 
and sleptons, and the potential is given by [l2| 



m2|$|2, (|$|« M s ), 

(1) 



(log§£) , (|0>|»M S ), 



in the gauge- mediated SUSY breaking 1 . Here ~ 
O(TeV) is a soft breaking mass, and Mp and Ms are 
related to the F and A components of a gauge-singlet 
chiral multiplet S in the messenger, respectively, as 



Mp = 



(M 



M S = (S), 



(2) 



where g stands for the gauge coupling in the standard 
model in general. Mp ranges as 



a 1 / 2 / 

10 3 GeV < Mp < ^— y / m 3/2 M P , 



(3) 



1 Mp and are related to parameters used in Ref. [Tl[ as Mp 

2~ 1 < /4 (m s M m / g) 1 / 2 and = V2m s , respectively. 
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where TO3/2 is the gravitino mass, and Mp(ss 2.4 x 10 18 
GeV) is the reduced Planck scale. 

Q balls form from the condensate of the flat direction 
during its helical motion after it started rotation when 
H ~ Mp/<f> osc , where (j> osc is the amplitude at the onset 
of rotation. The charge of the created Q ball is estimated 
as 



smaller than Q cr , where 

Qcr = 47T 4 C 4 



Mp 



TiNLSP 



(11) 



III. Q-BALL DECAY 



(4) 



The decay rate of the Q ball has an upper bound as 
111 4 . 



where (3 ~ 6 x 10~ 4 for a circular orbit (e = 1), while 
f3 ~ 6 x 10 -5 for an oblate case (e < 0.1). Here e denotes 
the ellipticity of the field orbit. The charge Q is in fact 
the $-number, and relates to the baryon number of the 
Q ball as 



B = bQ. 



Here we assign the baryon number b to $-particle. For 
example, b = 1/3 for the udd direction. The mass, the 
size, the rotation speed of the field, and the field ampli- 
tude at the center of the Q ball are expressed in terms of 
the charge Q respectively as 



Mr 



(M F Q 



3/4 



V2Tr(M F Q- 1/i , 



(6) 
(7) 

(8) 
(9) 



where £ is O(l) parameter determined by the fit to nu- 
merical calculation [TTI. [l3j 2 . 

The Q ball can decay into those particles kinemati- 
cally allowed. This condition is written as ujq > Wdecay 
Since we are interested in the production of baryons, it 
is rephrased as ujq > bm jf, where ~ 1 GeV is the 
nucleon mass, or ujq > m q before the confinement, where 
m q is the quark mass. In this article, we adopt the bot- 
tom mass so that m q — rrib ~ 4.2 GeV 3 , although it 
depends on when the Q balls decay and which flat direc- 
tion it is, for example. This leads to Q < Qd with 



r Q <T 



(sat,d) 

Q 



1 

Q96tt 2 



(12) 



derived by considering the maximum outgoing flow of 
the fermions x, which has coupling £j nt = /0*x?/+h.c, 
in the limite of Rqloq — > 00, and is valid typically for 
(5) f4>Q ^ W Q- Here we put the index d, since the elementary 
process is the decay <j) — > X + 7 ?- This rate is applicable to 
the Q-ball decay into NLSPs, and it actually expresses 



as r 



(NLSP) 

Q 



(sat,d) 



H3- 



On the other hand, the elementary process of the Q- 
ball decay into quarks when ujq < ttinlsp takes place 
through heavy gluino (or higgsino) exchange as <f> + <f> — > 
q + q. Since quarks can possess twice larger momenta 
than in the decay case, ujq in the decay rate should be 
replaced by 2ujq, so that, for the scattering, we have (Tl| 



^(sat,s) 

Q 



1 U Q 
Q \2tt 2 



AnRj 



(13) 



Numerical calculation in Ref. [Tl[ provides the prefactor, 
and the decay rate into quarks is given by 



r 



(q) 
Q 



1.1 N q r 



(sat,s) 



(14) 



where N q is the possible degrees of freedom of quarks, 
and we set it as 3 x 3 x 2 = 18 in the following. This is 
indeed the Q-ball decay rate that is used for the decay 
into baryons. 

In the opposite case of f<f>Q <C ujq, the rate was esti- 
mated in Ref. 14] as 



Q 



(15) 



4tt 4 C 4 ( ^ 
nib 



(10) 



Similarly, the Q ball decays into NLSPs, once the con- 
dition ujq > ?nNLSP is met. This leads to the charge 



for Rquiq — > 00, whose 
derstood as follows 
the surface layer of the C 
the penetration length of 
Therefore the decay rate 
uct of one-particle decay 



physical meaning can be un- 
The interaction region is only 
) ball whose depth is given by 
produced particles ~ (/</>q) _1 - 
can be estimated as the prod- 
rate (~ / 2 wq), particle density 
(~ wq^q), and the effective volume (~ Air Rq{j '^q) x ) ■ 



[TJ. In 



2 The actual value is estimated as ( = 2 1 / 4 c/n, where 
4.8 log(m s /u>o) + 7.4 in terms of parameters used in Ref. [l^ 
our successful scenario, typically c RJ 14, so that £ rs 2.5. 

3 Results do not change even if we take the top quark mass, since 
this condition is less restrictive. See Sec. V. 



4 The decay rate was first derived in Ref. Il4l . It is different from 
the present result by factor 2, since they show the pair production 
rate in the theory of x = V- 
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However, it was pointed out in Ref. [ll],[l3[ that this is 
not the case for realistic Q balls in SUSY. In fact, the size 
of the Q ball is Rq ~ Wq 1 , so that the penetration length 
of the produced particles becomes larger than the Q-ball 
size: (/^q)" 1 > Wq 1 ~ Rq. Therefore, the effective 
volume should be the whole volume of the Q ball. Then 
the decay rate has the following dependence: 



f$0\ p(sat,d) 

uq 



(16) 



The Q-ball decay into gravitinos is then given by [ll| 



p(3/2) _ q g ^UQ<t>Q ^ " p(sat ,11 



(17) 



where (F) — v3w3/2-Mp is the SUSY breaking F term. 



IV. ABUNDANCES 

Now we have in our hand the revised decay rates into 
baryons (|14p and into gravitinos p7|) . Since the main 
decay channel is the decay into baryons, the temperature 
at the decay is given by 



- (- 



90 



V47T 2 iV D 



1/4 



o.,c ; ,v,iV /2 ^ 1/2 ^ Vli 



Mp 



10 7 GeVy \10"' 



18 J V 61-75, 

V2 / n \ -5/8 

21 I > ( 18 ) 



where iVn is the relativistic degrees of freedom at the 
decay time. 

Refining the estimates of the baryon and gravitino 
abundances together with NLSP abundance, which are 
produced by the Q-ball decay, we have the relations be- 
tween those number densities and that of the (f> field as 

1 



On the other hand, the branching into NLSPs is derived 

as 



^(NLSP) 



-Bnlsp 



r 



(q) 



(23) 



once it is allowed kinematically when Q < Q CI . If the 
charge Q is smaller than Q cr when the Q ball formed, 
one should set Q C r/Q = 1 in Eq. (l2"Tl) . 

The actual abundances are dependent on whether or 
not Q balls dominate the energy density of the universe 
before the decay. The baryon abundances for the Q-ball 
dominated (QD) and non-Q-ball dominated (NQD) cases 
are thus given by 



Y h 



rib 
s 



37b n b 
4 PQ 



3Tp n b 
~ ~PQ 



3Td ebrid 



osc 

9T D eb 



4 |cjq?10 



37r,h n b 

4 /Orad 

9 



8^2 



RH 

ebf3~ 



16ojq 
3Trh n b 

4 Pinf 



3/4 MfTrh n3/4 

Ml 



(QD) 



where Eci.tfTS"]) and 

PQ 



M Q n Q ~ Mq— ~ -WQUtj,, 



Q 3/ \ (NQD) 
(24) 

(25) 



are used in the QD case, while ~ m c ff0 2 sc , rrieff = 
yJV" ~ 2\/2Mp / <f> osc , 3H osc = m off , and Eq.® are used 
in the NQD case. On the other hand, the abundance of 
the gravitino dark matter is determined by the ratio with 
the baryon density as 



(26) 



P3/2 _ m 3/2 "3/2 ^ m 3/2 ^3/2 
p b mjv n b TO at £& 



Therefore, we can eliminate e& in Eq. (f24|) using Eq.([26|) 
to obtain the baryon abundance as, setting (3 = 6 x 10~ 5 , 



n b ~ e6ri0, 
"-3/2 - B 3 / 2 n^, 

D Qcr 
"NLSP — -£>NLSP— ?p"0, 



(19) 
(20) 

(21) 



where -B3/2 — /Iq is the branching ratio of the 

decay into gravitinos estimated as 



B 



3/2 



I 2.5/ V 18 



f Mf Vf m ^ Y 2 (22) 
V10 7 GeVy V100 MeV/ ' v ; 



V (QD) / t N 7/2 / „ v -1/2 f N \ -1/4 



10 



-10 



2.5 



9 

18 



and 



/ ^3/2 y 

VIOO MeV/ 



[ 10 7 GcV / I 10 22 



(NQD) 



10 



-10 



61.75 

-3/8 



l.ii I) 5 ( — 



xf W3/2 V 
VIOO MeV/ 

X I 10 7 GeV J 



18 



10 4 GeV y 

3/4 



10 2 



(27) 



(28) 
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for the Q _ ball dominance and non-dominance, respec- 
tively. These relations must be hold for simultaneous ex- 
planation of the baryon and gravitino dark matter abun- 
dances. 

The NLSP abundance produced by the Q-ball decay is 
estimated as 



(Q) 
Pnlsp 



m 3 / 2 Y 3 / 2 



Pnlsp 



5m nY, 



5m N Y b 



P3/2 

TONLSP -SnlspQc 



^NLSP "-NLSP 
"3/2 



m 3/2 B 3/2 Q 



0.20 GeV 



10- 10 

WNLSP 



/ ™NLSP \ 

V300 GeV/ 



W3/2 
100 MeV 

3 / Q 



10 



22 



(29) 



at the decay time. Here and hereafter we take ttinlsp — 
300 GeV, although we keep showing its dependence. It 
is shown for some fixed value of the charge Q in the dark 
green dotted-dashed lines in Fig. [1] 

The NLSP abundance is usually determined by 
the subsequent annihilations of NLSPs as "nlsp — 
H(Tn)/(<Jv) The annihilation cross sections for the 
bino, stau, and sneutrino NLSP are estimated as UM 



(av)r r - 
(av)z c 
They lead to 

(ann) 

Pnlsp 



— ^aSfgv- 



(30) 
(31) 
(32) 



2.0 x 10" 4 GeV 



10.75 



-1/2 



( 



(ann) 

Pnlsp 



WNLSP \ 3 / Tp 

V300 GeV/ \MeV / 



(33) 



3.9 x 10" b GeV 



(- 



(ann) 

Pnlsp 



™NLSP 

V300 GeV 
1.2 x 10~ 6 GeV 

"T-NLSP 



V 10.75 

( T D 



VMeV 



(34) 



-1/2 



10.75 

3 / T D 



("3011 GeV 7 \\UY 



(35) 



When the density of the NLSP produced from Q-ball 
decay is larger than the annihilation density, Pj^gp > 

PnlsP' tne NLSP density will settle down to the annihi- 
lation value in the end. The annihilation abundances are 
displayed for several different decay temperatures Td in 
blue solid lines in Fig. [TJ 
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FIG. 1: NLSP abundances produced by the Q-ball decay for 
Q = 10 26 , 10 28 , and 10 30 in dark green dotted-dashed lines. 
The abundances fixed by the annihilation is denoted by blue 
solid lines. Black dotted line is the abundance that the NLSP 
decay gives the right amount of the gravitino dark matter. 
Red dashed lines are the BBN constraints [l6l ]. 



The abundance is generally limited by the BBN bound 
for each species of NLSP [lg, shown in red dashed lines, 
or bounded by the NLSP abundance that accounts the 
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right amount of the gravitino dark matter by the NLSP 
decay, denoted by black dotted line. The latter limit is 
given by [1, [l6| 



Pnlsp 



«3/2 



1.5 x 1(T 6 GeV ( 
100 McV 



™NLSP 



V300 GeV/ 
1 f Y h 



10 



-10 



(36) 



As we see later, the Q-ball charge is < 10 26 , so that 
the NLSP abundance is fixed by the annihilation den- 
sity. Therefore, for larger decay temperature T^, which 
will be achieved in the successful scenario, the gravitinos 
produced by the decay of NLSPs can entirely be neg- 
ligible for the gravitino dark matter abundance. More- 
over, severe BBN constraints are evaded, and we can take 
larger gravitino mass m 3 / 2 for successful scenarios. 

V. CONSTRAINING MODEL PARAMETERS 

Our scenario simultaneously explains the right amount 
of the baryon asymmetry and the gravitino dark matter 
both from the Q-ball decay. Since Y\, ~ 10~ 10 , for suc- 
cessful scenario, we can obtain, from Eqs. (f27| and (l28l) . 
the relations that hold for the charge Q as 



Q ~ 4.7 x 10 



17 



18 



2.5 

m 3/2 \~ 8 / 3 / Mp 



100 MeV 



61.75 

28/3 



V 10 7 GeV 



(37) 



Q ~ 2.4 x 10" 8 ( 



18 



rr \ - 4 / 3 / i\ r \ -20/3 

J- RH \ / Mp ■ 



10 4 GeV J \ 10 7 GeV J 



100 MeV 

1 



, (38) 



for QD and NQD cases, respectively. The latter is depen- 
dent on the reheating temperature after inflation, Trh- 
It has un upper bound so that thermally produced grav- 
itinos do not overdominatc the universe, which results in 



T RH £ 3 x 10 6 GeV ( 



^3/2 

100 McV 



(39) 



Taking account of this upper bound in Eq. (|38p . we obtain 
the lower limit as 



Q ~ 1.2 x 10^ — 
1 2.5 



(Y w/3 (nX /3 ( m f 



1 s y . eV") 



-20/3 



(40) 



Notice that it is independent of m 3 / 2 . Eqs. (|37l) and (j40|) 
are shown in blue solid and red dashed lines, respectively, 
in Fig. H 

Together with the upper limit ((3|) on Mp, shown in 
magenta double-dotted dashed lines in Fig. [3J we have 
the allowed regions in the parameter space {Mp, Q). We 
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FIG. 2: Allowed region for m 3/2 = (a) 100 MeV and (b) 500 
MeV, which is surrounded by the Yi, line for QD case (blue 
solid line) and for NQD case with Trh being fixed (red dashed 
line), and by Mf limit (magenta double-dotted dashed line). 
Iso- Td contours are shown in dark green lines. We also show 
the kinematical limit (cyan dashed) and the free streaming 
constraint (green solid). Black long dashed line represents 
the trace of the upper limit of the allowed region when m 3 / 2 
changes. 



show in Fig. [2] two cases with different gravitino masses: 
(a) 100 MeV and (b) 500 MeV. 

The boundary of the Q-ball domination/non- 
domination is derived from the relation 



1 = 89 



-1/2 



M F 



3/2 



10 4 GeV J V 10 7 GeV J 



1/4 

61.75, 

9/8 

22 / > ( 41 ) 



10 



where we use Eqs. ([25)) . (fT9[) . and the lower equation in 
Eq . (|2"4"]1 . This leads to the following boundary of the Q- 
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ball dominance on (Mp,Q) parameter space, given by 



2.a 



18 



-8/9 



Mp 

in 1 G.-Y ! V 10 7 GeV ) 



61.75 

\ -4/3 



-2/9 



(42) 



above (below) which corresponds to the (N)QD case, 
shown in orange dotted-dashed lines in Fig. [5J 

As Trh decreases to lower values, the baryon abun- 
dance line in NQD case sweeps the hatched region. We 
show two cases of Trh/Tr^ 1 "^ = 10~ 2 and 10~ 4 in 
Fig. 2(b). 

The allowed region for m 3 / 2 = 100 MeV shown in 
Fig. [5{a) applies to stau and sneutrino NLSPs, while it 
does not to the bino, as can be seen from Fig. [1] that 
the decay temperature should be > 10 GeV for the bino 
NLSP in order to evade the BBN constraint, while > a 
few MeV is possible for stau and sneutrino NLSP. Simi- 
larly, the allowed region for m 3 / 2 = 500 MeV is applicable 
to stau and sneutrino cases, while the bino case is allowed 
only for Tp> > 16 GeV (below the dotted line). (See also 
Fig. QJa) where we show Tp> — 16 GeV in dotted line.) 

Black long dashed lines show the trace of the upper 
bound of allowed regions when m 3 / 2 is changed. There- 
fore, the whole allowed region is below this long dashed 
line and above red short dashed line, and the lower limit 
of the gravitino mass is obtained as m 3 / 2 ^ 42 MeV. 
On the other hand, one can take larger gravitino mass 
of O(GeV) in the stau case, and even much larger in 
the sneutrino case, as can be seen from Fig. [T] It is al- 
lowed typically for Td <L 10 GeV, which corresponds to 
the lower part in {Mp, Q) plane. It will eventually con- 
strained by the fact that such small charge Q balls may 
evaporate in the thermal bath, or the Q ball cannot be 
the gauge-mediation type, where the logarithmic term 
can no longer dominate over the mass term (See Eq.([T])). 

Now let us check that the allowed regions obtained 
above are consistent with more general conditions, which 
indeed turns out that they are less restrictive. We show 
the decay temperature Tp> in dark green lines, the kine- 
matical bound ujq > m q in cyan lines, and the free 
streaming limit in green lines in Fig. [5J 

First we see the constraint on the decay temperature 
Td ■ Since the decay of the Q ball should take place before 
the BBN, the decay temperature must be larger than 
~MeV. The Td contour is given by 



5.3 x 10 



2d 



4/5 



2.5, 

rp \ -8/5 

MeV J 



V, 
18 



4/5 



Mp 



( iVp 
V 10.75 

4/5 



-2/5 



10 7 GeV 



(43) 



Next we consider the kinematical bound. The Q ball 
is kincmatically allowed to decay into quarks. This is 
achieved for Q < Q D (Sec Eq. fTU)) ): 



Q < 5.0 x 10^ y -V 



cW m f y 

2.5/ \10 7 GeV J 



(44) 



In addition, since the gravitino mass is relatively small, 
the produced gravitinos could free-stream to affect large 
scale structures. In order to avoid such a case, we impose 
that the free streaming length should not exceed ~ Mpc 
[18], which results in the bound on the present day free 
streaming velocity as vq < 3.7 x 10 -7 . vq is estimated as 



1 T / N 

Vq — ~ 

2 m 3 /2 T D \N D 



1/3 



(45) 



where To and Nq are respectively the temperature and 
the relativistic degrees of freedom at present. Therefore, 
the constraint is obtained as 



Q < 3.1 x 10 



2!) 



C V 4/3 fN r ^ 4/3 



2.5 

m 3 / 2 



V, 



D 



/ m 3/2 \ 

V100 MeV/ 



s/3 / Mp \ 
^ 10 7 GeV/ 



10.75 

-4/3 



2/9 



(46) 



VI. CONCLUSIONS 

We have reinvestigated the scenario that the amount of 
the baryons and the gravitino dark matter is naturally ex- 
plained by the decay of the Q balls in the gauge-mediated 
SUSY breaking. Using the more correct decay rates into 
gravitinos and baryons derived recently in Ref . [ill ] , and 
also taking into account the Q-ball domination and/or 
the free streaming condition, we have found that the sce- 
nario actually works. Although the branching into grav- 
itinos is smaller than the previous estimate, which leads 
to the more relative production of NLSPs by the decay 
of the Q balls, the actual abundance of the NLSP is sup- 
pressed by the annihilation effects. 

We have found that the scenario works typically for 
m 3/2 = 40 MeV - 5 GeV, Mp = 3 x 10 7 - 3 x 10 8 GeV, 
and Q — 10 18 — 10 24 , although there are some differences 
among the NLSP species. Therefore, it is natural to have 
this scenario of simultaneous direct creation of the baryon 
asymmetry and the gravitino dark matter from the Q-ball 
decay. 
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Appendix A: Indirect production of the gravitino 
dark matter via NLPs 



We briefly summarize the case of the indirect produc- 
tion of the gravitino dark matter from the Q-ball decay. 
The Q ball first decays into NLSPs, which eventually de- 
cay into the gravitinos that become the dark matter of 
the universe. In actual scenario, there are two cases. The 
one is that the NLSPs annihilate to reduce their abun- 
dance to the annihilation density fEqs. (|33|) — (f35l)L This 
is essentially the same case considered in Ref.(l(|. The 
other is the case that the amount of the NLSPs produced 
by the Q-ball decay is small enough that they do not an- 
nihilate to reduce the amount. In either case, the NLSP 
abundance should coincide with Eq. ()36[) . 

In the former case, where the annihilation is effective, 
it can be readily read off from Fig. [T] for the allowed 
range of the gravitino mass TO3/2 (or equivalenlty, the 
decay temperature of the Q ball, Tb): 



m 3/2 = 2.8 - 20 MeV, 
(T D = 3.7 - 26 MeV) 

m 3/2 = 39 - 500 MeV, 
(T D = 1 - 13 MeV) 

= (3 - 5) x 100 MeV. 
(T D = 1.6 - 3.9 MeV) 



(bino) 
(stau) 
(sneutrino) 



(Al) 



Notice that the lower bound in the bino NLSP is deter- 
mined by Tq^ > TnlsPi where the NLSP decay rate is 
estimated as 



NLSP 



NLSP 



487rm 2 /2 M 2 • 



(A2) 



As for the baryon abundance, it is determined by Eq. (j2~4")) 
with appropriate value of e (and Trh). 

On the other hand, the latter case appears in rela- 
tively narrow range of To- The argument is similar to 
the direct production discussed in Sec. IV and V. Since 
the gravitino abundance has straightforward relation to 
the baryon abundance through the number density of the 
NLSPs as 

-Bnlsp 



5 _ P3/2 _ 1TI3/2 n 3 / 2 ^ m 3/2 

Pb mpf rib TUN Q 
the baryon abundance is estimated as, from Eq. ([2~4 



(A3) 



Y 



(QD,indir) 



10 



-10 



9-") 10" ( X 



Nr 



10.75 



r /4 ( 



v 7/2 

2.5 J 

WlNLSP 



18 



V 



10 



Mf \ 

7 Gcvy 



300 GeV/ 

7/2 



Q_ 

22 



10 



-1/2 



W3/2 

20 MeV 

-11/8 



(A4) 



when the Q balls dominate the energy density of the uni- 
verse, and, for the non Q-ball domination case, 



Y 



(NQD,indir,s) 
IF 15 



0.42 



2.5 



18 



m NLSP 



x(- 

V 300 GeV J V 20 MeV 

( Trh \ / Mf \ 
X ^ 10 4 GeV ) \ 10 7 GeV J 

for ujq < mNLSp(= 300 GeV), and 

r(NQD,indir,b) 



10 2 



-1/4 



,(A5) 



Y; 



10 



-10 



0.23 



§ ( 



m 3 / 2 
20 MeV 



3/4 



Tkh \( Mf \( Q V 
10 4 GeV J \ 10 7 GeV J \ 10 22 J ' { ' 



for ujq > mjMLSP- We therefore obtain the relation be- 
tween Q and Mp from these estimates with Yb = 10 -10 . 
They are shown as an example of the stau and sneutrino 
NLSP for m 3/2 = 500 MeV in Fig. El Thick solid lines de- 
note the case when the annihilation fixes the NLSP den- 
sity, while the Q-ball decay determines the NLSP density 
for the hatched regions. 



28 



27 



10 26 
10 25 
O 10 24 
10 23 
10 22 



21 



20 













--J. \ NVeV 








— r^TTN^ 






















4 











10 J 



10° 10' 
M F [GeV] 



10° 



FIG. 3: Allowed region for m 3/2 = 200 MeV. The NLSP (or 
equivalently the gravitino dark matter) density is fixed by the 
annihilation (thick solid lines) or determined by the Q-ball 
decay (hatched regions). Red and blue regions correspond to 
the stau and sneutrino NLSPs, respectively. 



Finally, we comment on the free streaming constraint 
in the case of the indirect production the gravitino dark 
matter. Since the decay is determined by the NLSP de- 
cay rate (| A2|) . the present-day free streaming velocity is 
estimated as 



1 tonlsp 



To 



2 m 3 /2 Tnlsp 



N D 



1/3 



(A7) 



where Tnlsp = -T NL pg. This turns out the constraint on 
the NLPS mass as m NLSP > 270 GeV. 
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